An exactly solvable model of the persistent currents in presence 
of hard-core interaction in a one dimensional mesoscopic ring 
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A simplest model of the persistent current in presence of many-body interaction in a mesoscopic 
ring is presented. The Bethe ansatz approach is used to find exact solutions. Both cases in the 
absence of impurity and in the presence of impurity are solved and compared. It is indicated that 
both impurity and hard-core interaction do not affect the result of persistent current explicitly if 
there is a large number of electrons in the ring. 
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ON . It is well known that quantum coherence plays a central role in mesoscopic physics. The persistent current on 
mesoscopic rings threaded by a magnetic flux is a particular sensitive probe of such coherence. The existence of 
persistent currents on mesoscopic rings was predicted early in [1-3], but it had not absorbed much attention for a 
long time due to the technological reasons. Owning to the advances in nanotcchnology, more research interests were 
involved in both theoretical [4,5] and experimental [5-7] aspects. Among these theoretical researches, || studied 
the model that the electrons move freely in an ideal ring. The effect of a periodic potential was considered in the 
framework of hopping model Q . The influence of impurity was calculated mumerically [ |i"cfl and analytically in the 
framework of the grand canonical ensemble [ pd| The consideration of a screened Coulomb interaction was discussed 
[12 1 via the Hartree-Fock approximation, and the long-range part of the Coulomb interaction was taken into account 
[13 1 by means of the energy-density functional. There are also other discussions, such as a combination of impurity- 
t^. . and Coulomb- perturbation [ pd|JT^ ] , hopping models with a short-range two-body interaction |l(|[l7j , and the model in 
terms of the Luttinger liquid |iqor that from the Bogliubov-de Gennes equation |l]| etc.. However, as far as we were 
aware, there is not a thorough discussion about this subject by an exactly solvable many-body model. Present letter 
is an attempt toward this destination. As a result of our discussion, it is indispensable that the wave-function must be 
piecewisely continuous functions if one consider N electrons in a one dimensional ring. As long as the circumference 
of the ring is supposed to be the mesoscopic scale such that the inelastic scattering is considered not to happen, it is 
I i applicable to use the Bethe ansatz approach to solve the present problem. In the following, we discuss the problems 
of N particles and that of N particles with a single impurity on a mesoscopic ring in the presence of magnetic flux 
respectively. The results of both cases are compared. 

The model Hamiltonian for the system is supposed to be a A-body problem with delta-function interactions on a 
ring. For clarity and convenience, we convert the nature of the ring to be a system on a line with periodicity of L. 
^ ' The Hamiltonian reads 

^ '. tt 9 .2ir (j) n2 c, rs 

where u stands for the coupling constant; n, </>o and <f> stand for integers, the flux quantum hc/e and the total magnetic 
flux threading the ring, respectively. 

It is easy to show that if ip' is an eigenfunction of the Schrodinger equation in the absence of magnetic flux, 

then tp(xx,X2, ■ • •) = exp [ — '- — {x\ + x% H ) ) ip'(x\,X2, ■ • ■) satisfies the Schrodinger equation in the presence of 

magnetic flux. Moreover, if ip 1 satisfies the periodic condition then ip does not satisfy periodic condition but an 
additional torque (non-integrable phase) exp(i27T(/)/0o) comes about. However the wave-function of the real system, 
being on a ring, is required to be single valued. This implies that the periodic condition must be imposed upon the 
ip i.e. . ip(xi, ■ ■ • , Xi + L, • ■ •) = ip(xx, ■ ■ ■ , Xi, ■ ■_). Apparently, it is guaranteed as long as the wave-function ip' is 
assumed to satisfy a twisted-periodic condition J2y,|l[ instead of the usual periodic one. Therefore, the calculation of 
the energy spectrum of a Schrodinger equation in the presence of a magnetic flux with periodic condition is equivalent 
to that in the absence of the flux with, whereas, a twisted-periodic condition, i.e. 

ip(xi, ■ ■ ■ ,Xi + L, ■ ■ •) = exp(i2Tr(p/4>o)ip{xi, ••■,£*,•■•) (2) 

Because of the above mentioned equivalence, we only needs to solve the eigenequation of Hamiltonian (Q) by taking 
the flux as null while taking account of the twisted periodic condition (0). We emphasize that the characteristic 
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length L of the circumference is within the mesoscopic scale such that the inelastic scattering does not happen. This 
suggests the validity of Bethe ansatz approach for present problem because the ansatz embodies the property of non- 
diffractive (i.e. elastic in one-dimension) scattering. Apart from the consideration of the twisted periodic condition, 
the whole course of solving the problem is well known [^2],^3| . So we do not exhibit the details but only give the main 
results that are tightly related to our following discussion. The energy spectrum of Bethe ansatz solution is given by 
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E = ^— kj, where {kj} are roots of the following transcendental equations 



kj = (I 



(3) 



The {Ij} play the role of quantum numbers, and take integer value or half- integer value accordingly as N — is odd or 
even. As the Bethe ansatz wave function becomes null as long as two fc's take the same value, the system has fermionic 
property even for N-bosons. Thus the quantum numbers of ground state are {— -^j-^, — ^5-^- + 1, • • • , ^-t^ — 1, ^p-} 

If let kj = kj — -y- ("t - )j (||) becomes the usual secular equation for {kj} 
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in which the magnetic flux is not involved (and it agrees with the result of p2fl). So the roots of (Q) are independent 
of the flux (j>, which is useful in deriving the following (^|). In terms of these roots we can write out the energy of the 
system as follows: 



E = 
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It is different from the result of the free electron model where kj is an integer instead of the roots of the secular 
equation (|j). The knowledge of the full energy spectrum enables us to calculate the persistent current of the system 

dE 

by using the formula in |Q. For a single ring at zero temperature, the persistent current is calculated by I ((f) — ~ C ~J^' 
explicitly, 
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In deriving the second line of (Q) has been used. As a result, n := jj 53j=i Ij ^ s an integer regardless of TV being 
even or odd. For <fi = 0, the ground state corresponds to n = 0, and the first excited state corresponds to n = ±1. 
However, for <fi = ±</>o, the ground state corresponds to n = =Fl. Similarly, one can find that the current I(4>) is a 
sawtooth- type function with periodicity of <f>a. 

Now let us turn to the discussion of the situation in the presence of a single impurity, the Hamiltonian is given by 
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where x' stands for the coordinate of the impurity and Xi for that of the ?th particle. It should be addressed again 
that the energy spectrum for the given Hamiltonian is equivalent to that of the same Hamiltonian in the absence of 
magnetic flux but with a twisted-periodic condition instead of the usual one. Similar to our previous discussion in 
p4|, we take a scalar transformation x' — > xq = \J m'/mx'. With this transformation, one can write the kinetic part 
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of both electrons and the impurity to a standard Laplace operator in an (N + 1) dimensional euclidean space with 
cartesian coordinates. At the same time we would obtain the boundary condition arising from the delta-function terms 
conveniently. By using Gaussian integral theorem, one can derive the correct discontinuity relations for the derivatives 
of wave-function arising from the delta-function in Hamiltonian (fjj). The total energy is E = |^ S^Li + ~^>^ an d 



2m' 



2m ^j=l 

the corresponding wave-function is obtained to be ip(yo> j/i, • • • , Un) — e lKyo (p(yi, • • • , 2/at), where K := X^=i kj + X 
and yi{i — 0, 1, 2, • • • , N) are the coordinates of particles in the frame of reference of impurity. The twisted periodic 
boundary condition leads to K = n'^f) + (N + l)^(^). The Bethe ansatz wave-function (p is a piecewisely 
continuous function defined on separate regions which are specified by permutation group |24| Sjv+i, namely 
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where ok stands for the image of a given k := (ki, &2j • • • ; fejv) by a mapping of permutation a G S^r and the coefficients 
A(<7, r) are functionals on Sn®Sn+i. It is worthwhile to mention that the summation in (§h runs over the permutation 
group SjV) whereas the various regions on which the wave- function is defined are specified by S/v+i- This differs from 
the Bethe- Yang ansatz p2]| . 

All the coefficients A's in (||) are determined up to an overall scalar factor. The coefficients among different regions 
are related by the following relations: 



A(<j,<TiT) = rjA(aia,r), fori = 1,2, ■ • • , iV — 1. 
in which 77 = 1 (or — 1) for totally symmetric (or totally anti-symmetric) state, and 

v - i[(ak)i - A//i] 



A(a, o n t) 



v + i[(<rk)i - X/n] 

where /i = m'/m; <7n G Sn+i but G" Sn- The relation of the coefficients on any one of the regions is given by 

M x M-i[(qfc) t - (o%fi]g , , 

A\<Ji<T, T) = — — Tj. 
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If y = (yi , • • ■ , • • ■ , j/at) is a point on the region specified by r G Sjy+i, as a result of periodicity, the coordinate y' = 
(j/i, • • • , yi+L, • ■ • , yjv) must correspond to a point on another region specified by 7T G Sjv+i, where 7 = ctatA and A := 
0jV-i • • • G <Sjv- Then the requirement of the twisted-periodic condition yields ipy T (y') = exp(i2ir(f)/(f)o)ip T (y). 
After writing out this relation in terms of the Bethe ansatz solution (|^) , we find that the twisted periodic condition 
is guaranteed as long as the relation A{a^T)e lf - crk ^ lL = exp(i27r<^/0o) j 4(°'j T ) is imposed. By making use of (|To|), ( |Tl| ) 
and (^) repeatedly, we obtain 



v — i[(ak)i — X/ fi] -q- u — irj[{ak)i — (crfc)j 



v + i[{ak)i - X/fj] u + i[{ak)i - (ak)j] 
Taking the logarithm of (|l2|), we have a system of coupled transcendental equations, 
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where Ij takes integer values or half integer values respectively for N being an even or odd number. (|13|) is the secular 
equation for the spectrum and the {Ij} play the role of quantum numbers. From (^) we learned that an eigenstate 
related to {hi, &2, • • • , kx} has the same eigenenergy as those related to {A, fe, • • • , /cat}, {fei, A, • • • , /cat} etc. if fi = 1. 
However, the N + 1-fold degeneracy is broken by the replacement of one of the particle with an impurity. 

~{^r) an< ^ ^' = ^ ~ ^^(t")' ( B ) becomes the same secular equation as in Q: 
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in which the magnetic flux is not involved. Thus the roots of the equation are independent of the flux (j>. In terms of 
these roots we can write out the energy as follows 
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The persistent current is formally calculated as 
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which is just the same result of (g) for N + 1 electrons. This result is of no surprise. As an equilibrium quantity, the 
persistent current at zero temperature is obtained by taking summation of the currents over all states with energies 
below the Fermi level. Due to the hard-core interaction, The Fermi statistics confers upon those states. The impurity 
plays a similar role as an electron in occupying those lowest-lying single-particle levels with the corresponding quantum 
number A', which is dependent on the mass of impurity. For the originally occupied state with k'j for the jth electron, 
the quantum number would change into what associated with impurity occupying the state. The current for any 
occupation of N electrons and an impurity in Fermi sea and at the Fermi surface is the same. However, electron 
beneath the Fermi surface is raised to the Fermi surface or becomes excited state, this increases the result in (||) 
by + ^). We would emphasize that the effect of impurity scattering have not considered here so that we 

have not dealt with the effect of localization in the amplitude of persistent current which is given by that of impurity 
slowing down the moving electron. 

Obviously, the assumption of the impurity having the same charge as the electrons might be too special. If the 
impurity has a different electric charge q, one can make the same discussion as we did in the above. In the case, the 
twisted-periodic condition reads 



ip{xi, Xi + L, xn,x') = exp(i2iT(j)/(})o)tp(xi, x,, xn, x' 
ip(xi, xn, x' + L) = exTp(i2m>4>/(j)o)'ip(xi, Xi, ...,xn, x 



where v — q/e. The formula for the persistent current becomes 
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Clearly, the variance of electric charge of impurity from that of electrons affects the persistent current only when the 
number of electrons is not very large. Finally it needs to indicate that the formula ([if]) and ([l^) are not valid for the 
case of /i = co. As far as this case, the problem will become iV electrons moving in a ring with a narrow ((5-function 
type) junction because the model Hamiltonian do not have translational invariance. 

In summary, we have proposed an exactly solvable model in which the interactions of both electron-electron and 
electron-impurity are taken into account. As a result, the persistent current in a one dimensional mesoscopic ring 
are exactly the same for either free electrons or electrons with hard-core interaction, though the energy spectrums 
of them might be quite different. On the other hand, from the above discussion we found that a single impurity do 
not affect the result of persistent current at least for the case of hard-core interaction if there is a large number of 
electrons in the ring. 
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